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ABSTRACT 

Context. The injection of hydrogen into the convection shell powered by helium burning during the core helium flash is commonly 
encountered during the evolution of metal-free and extremely metal-poor low-mass stars. Recent multidimensional hydrodynamic 
simulations indicate that the hydrogen injection may also occur in more metal-rich stars due to turbulent entrainment which accelerates 
the growth of the shell convection zone and increases its size. However, ID stellar models cast doubts that helium-flash hydrogen 
mixing does occur as it requires the crossing of an entropy barrier at the helium-hydrogen interface. 

Aims. With specifically designed multidimensional hydrodynamic simulations, we aim to prove that an entropy barrier is no obstacle 
for the growth of the helium-burning shell convection zone in the helium core of a metal-rich Pop I star, i.e., convection can penetrate 
into the hydrogen-rich layers for these stars, too. We further study whether this is also possible in one-dimensional stellar evolutionary 
calculations 

Methods. We artificially shift the hydrogen-rich layer closer to the outer edge of the helium-burning shell convection zone in a Pop I 
star with a mass of 1.25 M Q , and simulate the subsequent evolution in 2D and 3D, respectively. We also perform stellar evolutionary 
calculations of the core helium flash in metal-rich stars implementing turbulent entrainment by means of a simple prescription. 
These simulations were performed with the Eulerian hydrodynamics code HERAKLES and the stellar evolution code STAREVOL, 
respectively. 

Results. Our hydrodynamical simulations show that the helium-burning shell convection zone in the helium core moves across the 
entropy barrier and reaches the hydrogen-rich layers. This leads to mixing of protons into the hotter layers of the core and to a rapid 
increase of the nuclear energy production at the upper edge of the helium-burning convection shell - the hydrogen injection flash. 
As a result a second convection zone appears in the hydrogen-rich layers. Contrary to ID models, the entropy barrier separating 
the two convective shells from each other is largely permeable to chemical transport when allowing for multidimensional flow, and 
consequently, hydrogen is continuously mixed deep into the helium core. We find it difficult to achieve such a behavior in one- 
dimensional stellar evolutionary calculations. 
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1. Introduction 

The hydrogen injection flash is believed to commence during the 
core helium flash of metal-free and extremely metal-poor (metal- 
licity Z < 10~ 4 ) low-mass stars when their helium core shell con- 
vection zone penetrates into the hydrogen-rich layers (Fujimoto 
et al. 1990; Hollowell et al. 1990; Schlattl et al. 2001; Cassisi 
et al. 2003; Weiss et al. 2004; Campbell & Lattanzio 2008). Such 
an injection of protons is not expected for more metal-rich low- 
mass stars because the ignition of helium occurs further away 
from the hydrogen-rich layer and the entropy barrier between 
the helium and hydrogen layers is higher. 

Multidimensional hydrodynamic simulations demonstrate 
that entropy barriers in stars are not as impenetrable as they 
may seem to be at first sight. The mixing beyond conventional 
Schwarzschild or Ledoux convection boundaries and the growth 
of convection zones on dynamic timescales due to turbulent en- 
trainment are proves of it (Meakin & Arnett 2007). Recent hy- 
drodynamic simulations of shell convection during the core he- 
lium flash also suggest that convection could reach the outer 



hydrogen-rich layers of the star within a week at around the 
peak of the core helium flash even in metal-rich stars (Mocak 
et al. 2009, 2010). However in these models, the hydrogen- 
rich layers are still far away from the helium-burning shell con- 
vection zone preventing the simulation of the expected pro- 
ton injection process within an affordable amount of comput- 
ing time. Alternatively, one may study the evolution during the 
core helium flash up to a possible hydrogen injection by one- 
dimensional stellar evolutionary calculations which include the 
effects of turbulent entrainment by means of a simple entrain- 
ment law. As such simulations are still missing, we implemented 
in our ID stellar evolution code a simple entrainment law based 
on the ideas of Meakin & Arnett (2007), and performed stellar 
evolutionary calculations of the core helium flash in low-mass 



1 We find that the growth rate of the shell convection due to turbulent 
entrainment is of the order of a few meters per second for our hydrody- 
namic models, the exact value depending on the stability of the convec- 
tion zone boundaries and the convective velocities near these bound- 
aries (Meakin & Arnett 2007). The latter in turn depend not only on 
dimensionality of the flow (2D, 3D), but also on the spatial resolution 
of the simulation. 
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Fig. 1. (a; left panel): Temperature distribution of the modified initial model (MOD, solid) and of the original model (M, dotted) as 
a function of radius r and of enclosed mass M, in solar units M , respectively. The dashed curve gives the entropy profile of model 
MOD. (b; right panel) Mass fraction of the main isotopes of model MOD as a function of radius. The hydrogen shell corresponds 
to the shaded region. 



Table 1. Some properties of the modified initial stellar model 
MOD: Total mass M, mass Mn e and radius Rn e of the helium 
core, total nuclear energy production by helium burning L He , and 
by hydrogen burning L P H and lJ w where the superscripts p and / 
denote that the respective numbers were obtained using a partial 
and full CNO network, respectively. 
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Pop I stars. The entrainment rates adopted for these calculations 
are taken from three-dimensional hydrodynamic simulations of 
Mocak et al. (2010). 

Hydrodynamic simulations of similar phases of stellar evolu- 
tion are currently also carried out by Herwig et al. (2006, 201 1). 
The results of their and our studies can have important implica- 
tions for the hydrogen injection not only during the core helium 
flash but also during the early- AGB phase {i.e., the "dual shell 
flashes" described in Campbell & Lattanzio (2008)) or in more 
advanced nuclear burning stages. 

The paper is organized as follows. We introduce the numer- 
ical tools and the stellar input model used for our simulations 
in the next section. The hydrodynamic simulations of the hy- 
drogen injection phase and the stellar evolutionary calculations 
mimicking turbulent entrainment during the core helium flash 
are presented in Sect. 3, and Sect. 4, respectively. A summary of 
our findings is given in Sect. 5. 

2. Physical Input and Codes 

For the multidimensional hydrodynamic simulations we used 
the HERAKLES code (Mocak et al. 2008) derived from the 
PROMETHEUS code (Fryxell et al. 1991; Miiller et al. 1991). 
In brief, HERAKLES solves the Euler equations coupled with 
source terms describing self-gravity and nuclear burning. The 
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Fig. 2. Nuclear energy production rate as a function of radius r 
for the 3D model exphif.3d.f at t = s (dotted-black), and at 
t = 9800 s (solid-red), respectively. 



hydrodynamic equations are integrated with the PPM recon- 
struction scheme (Colella & Woodward 1984) and a Riemann 
solver for real gases according to Colella & Glaz (1984). The 
evolution of the chemical species is described by a set of addi- 
tional continuity equations (Plewa & Miiller 1999). Self-gravity 
is handled according to Miiller & Steimnetz (1995), gravita- 
tional potential being approximated by a ID Newtonian potential 
obtained from the spherically averaged mass distribution. The 
code integrates nuclear networks with the semi-implicit Bader- 
Deufelhard method (Bader & Deuflhard 1983; Press et al. 1992). 

As an initial model we used the helium core structure of a 
1 .25 M Q , Z — Z star during the core helium flash at its peak 
nuclear luminosity (Tab. 1). In the initial model, computed with 
the GARSTEC stellar evolution code (Weiss & Schlattl 2000, 
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Table 2. Some properties of the 2D and 3D hydrodynamic simulations based on model MOD: wedge size w; number of grid points 
and resolution in r, 9, andcp direction, respectively; estimated Reynolds number R c ; characteristic velocity v c of the convective flow 
in helium-burning layers; typical convective turnover timescale r conv at time at t ~ 40000 s for model 2d.p, and at t ~ 2500 s for 
models 2d.f and 3d.f, respectively; t max final evolutionary time. 
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2007), the hydrogen-rich layers are artificially shifted closer to 
the outer boundary of the helium-burning shell convection zone 
at radius around 9.8 x 10 8 cm. Originally the lower boundary 
of the hydrogen layer is located at a radius of 1.9 x 10 9 cm. 
The unmodified initial model (M) is described in more detail in 
Mocak et al. (2008, 2009). Compared to this model the resulting 
hydrostatic structure of the modified initial model MOD differs 
by a temperature sink at the boundary between the helium and 
hydrogen-rich layers (H-He boundary) at 9.5 x 10 8 cm < r < 
10 x 10 8 cm (Fig. 1). At this boundary the entropy increases al- 
most discontinuously by 20%, producing the so-called barrier. 
The temperature gradient in the helium-burning convection zone 
(4.7 X 10 8 cm < r < 9.2 x 10 8 cm) remains almost unaltered as 
well as the nuclear energy production. 

At the H-He boundary of model MOD around r ~ 9.8 X 
10 8 cm, we find a narrow peak in the nuclear energy produc- 
tion rate (Fig. 2) that is due to some left-over hydrogen (X ~ 
10" 4 ) when modifying the original initial model. In the lay- 
ers near the H-He boundary the temperature is relatively high 
(T ~ 5 x 10 7 K), and even if the hydrogen mass fraction is small, 
the nuclear energy production is significant. Although, this fea- 
ture is not seen in ID stellar evolution simulations of the hy- 
drogen injection flash (Hollowell et al. 1990), it may represent a 
somewhat more "realistic" situation, where the H-He interface is 
not a discontinuity but a moderately hot layer partially enriched 
with hydrogen. 

For the stellar evolutionary calculations of the core helium 
flash we implemented a simple turbulent entrainment law in the 
STAREVOL code (Siess 2006) and evolved a 1 M Q , Z = Z G 
star for different entrainment rates guided by our 3D hydrody- 
namic simulations. To achieve hydrogen injection even in metal- 
rich low-mass stars, we also tried the classical overshooting pre- 
scription based on the exponential decay of the convective veloc- 
ity field beyond the convection zone boundaries (Freytag et al. 
1996) with various efficiencies. 

STAREVOL is a ID lagrangian stellar evolution code able 
to compute the evolution of stars in the mass range (0.1 - 60 
M Q ) from pre-main sequence up to Neon ignition. The nuclear 
network includes 53 species coupled by 177 nuclear reactions 
(Siess & Arnould 2008), an accurate equation of state tested in 
the domain of very low mass stars (Siess et al. 2000), and all 
relevant neutrino contributions (for a detailed description, the 
reader is referred to Siess (2006). This code has been widely 
used to study pre-main sequence stars (Siess et al. 1997), rota- 
tional mixing (Palacios et al. 2006), AGB nucleosynthesis (Siess 
et al. 2002) and super AGB stars (Siess 2010). 

3. Hydrodynamic simulations 

We performed two 2D (models hifexp.2d.f and hifexp.2d.p) and 
one 3D (model hifexp.3d.f) hydrodynamic simulation using dif- 



ferent nuclear reaction networks. The last letter of the model 
name, / and p, refers to simulations performed with a full and 
partial CNO reaction network {i.e., neglecting reactions related 
to 13 N and 15 0). These runs start from the modified initial model 
MOD with the shifted hydrogen profile. Some general properties 
of these simulations are summarized in the Tab. 2, all of which 
were performed on an equidistant spherical grid ranging in ra- 
dius from 2 x 10 8 cm to 1.6 x 10 9 cm. The adopted grid resolu- 
tion is motivated by our earlier work (Mocak et al. 2008), where 
we found convergence of the results starting at a radial grid res- 
olution Ar ~ 3.7 x 10 6 cm and an angular resolution of ~ 1°. 
The boundary conditions in the radial direction are reflective, 
while periodic boundaries are imposed in the angular directions. 
For models hifexp.3d.f and hifexp.2d.f the reaction network in- 
cludes the following species: 'H, 3 He, 4 He, 12 C, 13 C, 13 N, 14 N, 
I5 N, 15 0, 16 0, I7 0, 20 Ne, 24 Mg, and 28 Si, coupled by reactions 
of the full CNO cycle and the triple-cr reaction. More details are 
given in the PhD thesis of Mocak (2009). 

For the simulation hifexp.2d.p we omitted the isotopes 13 N 
and ls O, the jS-decay reactions of the CNO cycle and the 
I2 C(p,y) 13 N reaction. Although this variant is unrealistic, it al- 
lows us to demonstrate that entropy barriers which are not sus- 
tained by sufficiently strong nuclear burning can be completely 
destroyed by the growing convection zone due to turbulent en- 
trainment. When considering all species of the full CNO cycle, 
the nuclear energy production at the H-He boundary is a factor of 
3 higher than in the partial network and leads to the formation of 
a shell convection zone in the hydrogen-rich layers (Fig. 3; mod- 
els hifexp.2d.f and hifexp.3d.f; see Sect. 3.1). No shell convec- 
tion zone appears in the hydrogen-rich layers of the 2D model 
hifexp.2d.p. 

3.1. Simulations with the full CNO reaction network 

Our hydrodynamic 3D and 2D models, hifexp.3d.f and hif- 
exp.2d.f (Fig. 3), are evolved with an accurate treatment of the 
CNO cycle on a numerical grid having the same number of zones 
in radial and angular directions. The latter allows us to compare 
the results of these simulations being biased only by the dimen- 
sionality of the simulation, but not by the grid resolution. 

In the 3D (2D) model, hifexp.3d.f (hifexp.2d.f), after 800 s 
(1000 s) convection fully develops in the layers located above 
the temperature maximum T max (r ~ 4.7 x 10 8 cm). As soon as 
the first convective plumes reach the H-He boundary, the dredge- 
down of hydrogen sets in. Up to ~ 1500 s the composition profile 
in both models is similar, with a steep increase of the hydrogen 
mass fraction X('H) exceeding 10~ 9 just below the H-He inter- 
face (Fig. 4). Subsequently and until t ~ 8000 s, X('H) remains 
almost constant in the helium-burning convection zone of the 
3D model, where 10" 10 < X('H) < 10" 8 . This is in contrast 
with the 2D simulation, hifexp.2d.f, where X('H) continues to 
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Fig. 3. Top row: snapshots showing contour plots of the modulus of the velocity |v| (in units of 10 6 cm s , left panel) and the 
hydrogen mass fraction (right panel) for the 3D model hifexp.3d.f at t~ 13300 s in the meridional plane (<p = 0°). Bottom row: same 
quantities but for the 2D model hifexp.2d.f. The dotted line separates helium- and hydrogen-rich layers (H-shell), and the dashed 
lines mark the boundaries of the computational domain. 



increase during this epoch. Between 8000 s and 12 000 s, the 
proton mass fraction increases in the helium-burning convection 
of both 2D and 3D models. 

At t ~ 12 000 s, the entropy barrier at the H-He interface 
in both models is steeper and almost twice as large as it was 
at the beginning of the simulations due to the enhanced nuclear 
burning. This increase in temperature and entropy at the inter- 
face slows down the injection of protons and the hydrogen mass 
fraction in the helium-rich layers decreases (Fig. 4). 

As protons are dredged down into the helium core, they are 
captured by 12 C via the 12 C(p, y) 13 N reaction, i.e., the hydro- 
gen abundance decreases sharply with depth. Nevertheless, some 
protons can reach layers with temperatures exceeding 10 8 K, 
where their lifetime against the 12 C(p, y) 13 N reaction is < 100 s 
(Caughlan & Fowler 1988). This is much shorter than their life- 
time outside the He-burning convection shell. As the convective 
plumes have velocities of 10 6 cm s (2 x 10 6 cm s _1 ) in the 3D 
(2D) model, hydrogen-rich gas is advected from the "colder" 
top of the He-burning convection shell to the aforementioned 
hot layers (located at r ~ 7.7 x 10 8 cm s -1 ) in less then 200 s 
(100 s). Our simulations show that this time is too short for the 
hydrogen to burn completely while being advected to layers with 
temperatures > 10 8 K(Fig.4). 

The differences between the X('H) distributions of the 2D 
and 3D model can be understood by realizing that the convective 
velocities in the 3D model are almost a factor of two smaller 
than in the 2D one. Hence, there is more time for hydrogen to 
burn during its transport in the 3D model than in the 2D one, 
i.e., hydrogen is able to reach the hot layers in larger amounts 
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Fig. 4. Temporal evolution of the (logarithm of the) hydro- 
gen mass fraction X('H) at three different radii [temperatures] 
around the top, middle, and bottom of the helium-burning con- 
vection shell, respectively: r x = 8.8 x 10 8 cm [T\ ~ 0.8 x 10 8 K] 
(dash-dotted), r 2 = 7.1 x 10 8 cm [T 2 ~ 1.2 x 10 8 K] (solid), and 
r 3 = 5.4 x 10 8 cm [T 3 ~ 1.5 X 10 8 K] (dotted). The upper panel 
(a) gives the results for the 3D model hifexp.3d.f, and the lower 
panel (b) for the 2D model hifexp2d.f, respectively. 
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Fig. 5. Radial distribution of (a) temperature and (b) entropy at 
t = Os (dotted) and t\ = 13 500 s for the 2D model hifexp.2d.f 
(dashed, blue), and the 3D model hifexp.3d.f (solid), respec- 
tively. The inserts show a zoom of the region around the H-He 
interface. 



in the latter case (Fig. 4). Besides higher convective velocities, 
the convective plumes arising at the H-He boundary also fill a 
larger volume fraction in the 2D model compared with the 3D 
one, where the plumes are smaller and narrower. 

The nuclear energy production at the H-He boundary and 
the impacting convective plumes change the initial properties of 
the H-He interface. At the start of the simulations, the bound- 
ary is characterized by a steep radial step in both the temper- 
ature and entropy profiles. However, when convection reaches 
the interface, the single step structure of the temperature and 
entropy profiles evolves into a two steps structure exhibiting a 
narrow entropy plateau between the steps (Fig. 5). This plateau 
reflects the formation of a secondary convection shell, located 
in the hydrogen-rich layers and powered by the CNO cycle. We 
find higher absolute values of both the entropy and temperature 
in the 2D model, hifexp.2d.f, compared to its 3D counterpart 
across the forming new convection shell. This can be explained 



by the stronger nuclear burning occurring in that model, which 
also leads to a faster growth of the hydrogen-burning convection 
sheU. 

Nuclear burning at the H-He interface steadily increases dur- 
ing both simulations, and eventually leads to the formation of a 
distinct second temperature peak (Fig. 5, Fig. 7). This behavior is 
well known from ID stellar evolutionary calculations of the core 
helium flash in extremely metal-poor stars, where the helium- 
burning convection shell splits after the injection of hydrogen. 
However, the splitting process seen in ID stellar calculations 
differs from that found in our multidimensional simulations. In 
ID simulations the splitting occurs due to the appearance of a 
temperature peak in the convective helium-burning shell at the 
location where hydrogen burning by the 12 C(p, y) 13 N reaction is 
proceeding the fastest. 

In our hydrodynamic models the appearance of a second 
temperature peak is preceded by a slow migration of the H-He 
boundary into the helium-burning convective shell, which forces 
it to retreat slightly (Fig. 7). Eventually, when the nuclear energy 
release and the temperature at the H-He boundary become suf- 
ficiently high, we observe the appearance of a new convection 
shell in the hydrogen-rich layers above the original one. Hence, 
what actually happens is not a splitting, but a "retreat" of the ini- 
tial helium-burning convection shell and the birth of a new one 
on top of it powered by burning of hydrogen. 

Our hydrodynamic simulations exhibit another major dis- 
tinct property: the absence of an impermeable radiative layer 
between the two convective shells. Mixing of nuclear species be- 
tween these layers is not prohibited, but occurs during the whole 
event (Fig. 6) with a decreasing efficiency towards the end of 
simulation. 

The evolution of the hydrodynamic models is likely rep- 
resentative of early phases during the dual core and He-shell 
flashes, and indicates that the entropy barrier does not inhibit hy- 
drogen dredge-down. However, this conclusion should be mod- 
erated by the fact that the entropy at the H-He boundary can in- 
crease by the rise of the nuclear energy production. Eventually, 
this entropy barrier may become so high that mixing ceases 
completely. What we have shown here is that, given our initial 
conditions, nothing seems to prevent H-injection, at least in the 
early stages. A similar behavior was also found by Herwig et al. 
(2011) in their 3D simulations of the hydrogen injection into 
the helium shell flash convection zone which helped them to ex- 
plain major features observed in spectra of the Sakurai's object 
(Duerbeck et al. 2000) (post-AGB object experiencing its last 
helium flash). Finally note that proton injection may have in- 
teresting consequences concerning the production of s-process 
elements as some 13 C can be produced which is a source of neu- 
trons (Campbell et al. 2010). 

We note, that during its final evolution the 2D model hif- 
exp.2d.f eventually experiences another unrealistic nuclear run- 
away when He ignites in the secondary convection shell (due to 
high temperature > 10 8 K at the H-He boundary) causing it to 
rapidly expand to the outer boundary of the computational do- 
main, when we stopped the calculation (Fig. 6 b). Concerning 
this unexpected behavior we have to emphasize that these hy- 
drodynamical simulations are strongly dependent on the initial 
conditions. In particular, we note that in the initial model MOD 
the layers at the H-He boundary are highly degenerate (degener- 
acy parameter iff » 10). 2 Such a large degeneracy favors a nu- 
clear runaway. Additionally, our multidimensional simulations 



2 Matter at the H-He boundary after proton injection (and some ear- 
lier expansion) in a ID population III model of 0.85 M computed by 
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Fig. 6. (a) Temporal evolution of the angle-averaged radial distribution of (the logarithm of) the hydrogen mass fraction for the 
3D model hifexp.3d.f, and (b) for the 2D model hifexp.2d.f, respectively. The dashed lines delineate the boundaries of the newly 
forming hydrogen-burning convection shell, and the dotted line marks the bottom of the helium-burning convection zone. 




Fig. 7. Radial distribution of (a) temperature, (b) hydrogen mass fraction, and (c) energy production rate at the H-He interface for 
the 3D model hifexp.3d.f at five different times: t = Os, t\ = 3230 s, t 2 = 6500 s, f 3 = 10400 s, and t 4 = 14 340 s. 



are likely to underestimate the star's expansion and subsequent 
cooling because of the reflective boundary conditions imposed in 
the radial direction. These numerical restrictions limit the cool- 
ing and favor the aforementioned helium runaway. 

3.2. Simulation with partial CNO network 

In the 2D simulation hifexp.2d.p we did neither consider 
the /3-decay reactions n N(J3~) l3 C and 15 O0S") 15 N, nor the 
12 C(p, y) 13 N reaction, which gives rise to a reduced energy pro- 
duction rate by roughly a factor of three at the H-He interface. 
It allowed us to demonstrates, that entropy barriers at convec- 
tion zone boundaries that are insufficiently sustained by nu- 
clear burning are penetrable and unable to prevent the growth 
of the convection zone due to turbulent entrainment on dynamic 
timescales. By neglecting proton capture on 12 C in the simu- 
lation, we obtain an incorrect hydrogen profile and a reduced 
nuclear energy production rate. That explains why we do not 



Simon Campbell (private communication) (Mocak 2009) has a degen- 
eracy parameter if/ as -2. 



observe here the formation of the secondary hydrogen-burning 
shell convection as in our hydrodynamic models with the full 
nuclear CNO reaction network (see previous subsection). 

In this model, helium-burning convection starts at t ~ 1000 s 
and extends throughout the whole initially convectively unstable 
shell as determined by the Schwarzschild criterion. The convec- 
tion shell is located between the radius of the temperature maxi- 
mum and the base of the hydrogen-rich layers. It is characterized 
by vortices, which are typical structures of 2D simulations due 
to the imposed axial symmetry. Their width is equal to the radial 
extent of the convection shell, i.e., r ~ 5 x 10 8 cm, which is also 
the width of the convection zone in the unmodified initial model 
(M). 

The shell convection reaches the hydrogen-rich layers 
shortly after the appearance of the first convective plumes. It 
initiates the dredge-down of protons deep into the helium core, 
and causes the outer boundary of the convection shell to cross 
the entropy barrier at the H-He interface (Fig. 8). The value of 
the entropy at this interface constantly decreases due to turbu- 
lent entrainment, which effectively shifts the entropy barrier to 
larger radii. 
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Fig. 8. Mean molecular weight (a) and entropy (b) at the outer 
edge of the convection zone of the 2D model simulated with 
the partial CNO network (hifexp.2d.p) at four different times: 
to = s (dotted), t\ = 46 600 s (dashed), t 2 = 94200 s (dash- 
dotted), and ? 3 = 142 000 s (solid red). 



The entrainment rate of the outer convection boundary is ini- 
tially relatively high, around 70 m s _1 , but after roughly 30 000 s 
it drops and eventually relaxes to an almost constant value of 
~ 4.5 ms 4 . The entropy barrier has a stabilizing effect be- 
cause it decreases the entrainment velocity. However, it does 
not prevent the advance of the convection zone into the H-rich 
layers and leads to an erosion of the entropy gradient. After 
170000 s the shell convection has moved ~ 8 x 10 7 cm into the 
H-rich buffer, meaning that almost 12% of the hydrogen-rich 
shell included into the simulation (limited by the outer radial 
grid boundary at 1.6 x 10 9 cm) is already entrained within only 
47 hours, or ~ 300 convective turnover timescales. 

In this simulation, protons are able to reach layers with tem- 
peratures as high as ~ 1.8 x 10 8 K. However, such a deep mixing 
of hydrogen into the helium burning layers is not realistic be- 
cause protons should have been captured by 12 C before. 



4. Stellar evolutionary calculations 

To mimic the effects of the turbulent entrainment in ID stellar 
evolutionary calculations, we experiment and implement a sim- 
ple entrainment law 

dM 

M E = ~^-u E (1) 
or 

where M is the mass of the shell at radius r and ue is the entrain- 
ment rate (Meakin & Arnett 2007). We simulated a 1 M , Z = Z G 
star up to the core helium flash. At first we use an entrainment 
rate of 6.7xl0 2 cm s -1 as found by Mocak et al. (2009) in their 
3D hydrodynamic simulations of the core helium flash. The en- 
trainment is switched on only shortly before the flash reaches its 
maximum intensity (maximum nuclear energy generated in the 
helium core). We find that the He-convection zone boundary is 
unable to cross the entropy barrier and penetrate into the hydro- 
gen shell. This result also holds when increasing the entrainment 
rate by one or even two orders of magnitude. 

This is partially due to the inherent time dependence of the 
mixing prescription. Indeed, for a given timestep Af, mixing will 
extend over a mass range given by AM e = M E At. Practically, we 
simply shift the convective boundary to the new mass location 
M e c d „f v + AM e . As the helium convection shell approaches the H- 
rich layer, the timestep which is constrained among other things 
by the rate of change of the nuclear luminosity, decreases thus 
preventing the injection of protons and a subsequent H-flash. 
Eventually, a (small) timestep is found that allows the simula- 
tion to go on and after the some time, the helium-shell convec- 
tion zone naturally quenches and H-injection is avoided. 

On the other hand, if we experiment with overshooting by 
using the exponential decay of the convection velocity in the ra- 
diative layers, mixing is treated as a diffusive process with the 
diffusion coefficient given by 

(2) 



nv exp(-|r - r edge \f(2fH P ) 



where Hp is the pressure scale height at the edge of the con- 
vective zone and / a free parameter. Compared to the previous 
prescription, this treatment shows two main differences : mixing 
is treated as a diffusive process so the penetration of the He- 
shell convection is "smoother" and the new convective bound- 
ary at t + Af does not explicitly depend on its location at time 
f and on the timestep. Our ID experiments indicate that at so- 
lar metallicity, a substantial overshooting is required to trigger 
an injection with / ~ 1. This injection then leads to the split- 
ting of the He-shell convection and after some structural read- 
justment to the pollution of the envelope. We report C and N 
enrichments by a factor of 2 or 3. Interestingly, we find that sub- 
sequently to this mixing event, the 12 C/ 13 C ratio has dropped to 
below ~ 10. Moreover, these He-core burning stars are carbon- 
rich (C/0~ 1.7) and show large 7 Li abundances (e{Li) « 3.5). 
Some of these features are found in J-stars but contrarily to ob- 
servations, our models show a too large 16 0/ 17 ratio (by a factor 
of two at least, and the C/N ratio is too low in our models). 

5. Summary 

We performed hydrodynamic simulations of shell convection in 
the core of a low mass star at the peak luminosity of the helium 
flash. In these computations, the hydrogen-rich layers were ar- 
tificially shifted down into the helium core. They allowed us to 
demonstrate a penetration of the helium-buming shell convec- 
tion into the hydrogen layers due to turbulent entrainment de- 
spite the presence of an entropy barrier between them. Resulting 
consequences are the following. 
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Hydrodynamic models computed with the full CNO reac- 
tion network leads to the formation of a secondary shell convec- 
tion zone in the hydrogen-rich layers. It is powered by hydrogen 
burning and is located above the already-present helium-burning 
shell convection. These simulations show (contrary to ID stel- 
lar calculations of similar evolutionary phases) that the mixing 
of protons down to the helium-burning convection shell is not 
completely inhibited, although it eventually slows down due to 
an increase of the nuclear energy production and the growth of 
the entropy barrier at the interface between the helium-burning 
and hydrogen-burning convection shell. 

In our hydrodynamic model computed with a reduced 
nuclear reaction network the energy production rate in the 
hydrogen-rich layers is smaller and we do not observe an ap- 
pearance of the hydrogen-burning convection shell. This simu- 
lation shows, that an entropy barrier not sustained sufficiently 
by nuclear burning does not act as an obstacle for the growing 
helium-burning convection shell due to turbulent entrainment. 
The entrainment rate of the outer rim of the helium-burning con- 
vection zone reaches a roughly constant non-vanishing value. 

Given our initial model of low-mass and metal-rich star, the 
aforementioned situations were difficult to reconstruct by one- 
dimensional stellar evolution code. 
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